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Variations in the D/H ratio of extended sightlines from FUSE 

observations^ 

Crist ina M. Oliveira^ and Guillaume Hebrard^ 
ABSTRACT 

We use new FUSE data to determine the column densities of interstellar D I, 
NI, 01, Fell, and along the HD 41161 and HD 53975 sightlines. Together 
with A^(HI) from the literature (derived from Copernicus and lUE data) we de- 
rive D/H, N/H, and 0/H ratios. These high column density sightlines have both 
log A'"(HI) > 21.00 and allow us to probe gas up to ~1300 pc. In particular these 
sightlines allow us to determine the gas phase D/H ratio in a hydrogen column 
density range, log A^(H) > 20.70, where the only five measurements available 
in the literature yield a weighted average of D/H = (0.86 ± 0.08) xlO~^. We 
find D/H = (2.14 ± °j^)xl0-5 along the HD 41161 sightline. This ratio is ~3(t 
higher than the weighted mean D/H ratio quoted above, for sightlines with log 
A^(H) > 20.70, while the D/H ratio for the HD 53975 line of sight, D/H = (1.02 
± Q2o)xlO~^, agrees within the la uncertainties. Our D/H measurement along 
the HD 41161 sightline presents the first evidence of variations of D/H at high 
A^(H). Our result seems to indicate that either the long sightlines that accord- 
ing to the deuterium depletion model are dominated by cold undisturbed gas 
where deuterium would be depleted onto carbonaceous grains occur at higher 
A^(H) than previously thought or that the clumping of low D/H values in the 
literature for the long sightlines has another explanation. Another possibility to 
explain the high D/H toward HD 41161 is local production of DI which cannot 
be ruled out. Both of the 0/H ratios derived here, (9.12 ± ?:^^)xlO-^ and (5.37 
± J'f4)xlO~^ (for HD 41161 and HD 53975, respectively) are higher than what 
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has been found by other authors in studies containing many sighthnes, although 
not inconsistent, in the case of HD 53975, with individual measurements along 
some sightlines. However, 0/H along HD 41161 is extremely high, and together 
with other high 0/H ratios based on FUSE data might indicate the presence 
of unknown systematic effects when deriving A^(0 1) along high column density 
sightlines. Finally, we derive (N/H)xl05 = 8.32 ± f^l and 5.07 ± \il and 
(D/O)xl02 = 2.29 ± '^oil and 1.91 ± for HD 41161 and HD 53975, respec- 
tively. In addition, the relatively high signal-to-noise ratio of the HD 41161 data 
allows us to place constraints on the /-values of some neutral chlorine transitions, 
present in the FUSE bandpass, for which only theoretical values are available. 

Subject headings: ISM: Abundances — ISM: Evolution — Ultraviolet: ISM — 
Stars: Individual (HD 41161, HD 53975) 

1. INTRODUCTION 

Deuterium is thought to be produced in appreciable amounts only in the Big Bang (Ep- 
stein ct al. 1976). Of all the light elements produced during Big Bang Nucleosynthesis (e.g., 
^Hc, ^Hc, ^Li), deuterium, easily destroyed in stellar interiors (astration), is the one whose 
relatively simple evolution makes it a particularly sensitive estimator of the baryonic density. 
Conversely, measurements of the cosmic microwave background, such as those performed 
with WMAP and other missions, can be used to infer the primordial abundance of deu- 
terium, (D/H)p. Analysis of the three year WMAP data (Spergel et al. 2006) yields (D/H)p 
in agreement with that derived by Kirkman et al. (2003), (D/H)p = (2.78 ± g:||)xlO-^ by 
taking the weighted average of five D/H mcasTirements toward QSOs in the range z = 2.079- 
3.572. Sembach et al. (2004) found D/H = (2.2 ± 0.7) xlO'^ for Complex C, a high-velocity 
low-metallicity cloud falling into our galaxy, which has presumably experienced more stellar 
processing than the gas seen toward QSOs. 

Precise measurements of the gas-phase D/H ratio in the ISM were first performed with 
the Co|?ernzctts satellite (e.g., Rogerson & York 1973) followed by the Hubble Space Telescope 
twenty years later {HST, e.g., Linsky et al. 1995), IMAPS (Interstellar Medium Absorption 
Profile Spectrograph, Jenkins ct al. 1999; Sonncborn ct al. 2000), and more recently the Far 
Ultraviolet Spectroscopic Explorer {FUSE, e.g., Oliveira et al. 2006, and references therein). 
Variations in the D/H ratio were observed early on (York 1983) and by different satellites 
(e.g., Jenkins et al. 1999). Lemoine et al. (1999) discuss some of the mechanisms that could 
be responsible for the observed variations, based on the D/H data available at the time. 
However only after the launch of FUSE in 1999, has a large enough sample of sightlines 
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with differing D/H ratios been accumulated, allowing the question of what causes the D/H 
variations to be addressed in a more systematic way. While the constancy of the D/H ratio 
in the Local Bubble (LB, corresponding to log A^(H) < 19.2, Sfeir et al. 1999) seems to be 
a well accepted fact ((D/H)lb = (1-56 ± 0.04) x 10^^ from the compilation by Wood et al. 
2004) two distinct explanations have emerged to explain the large scatter of gas-phase D/H 
measurements, in the Solar neighborhood outside of the LB. These measurements range 
from D/H = (0.50 ± 0.16) xlO-^ toward 9 Car derived by Allen et al. (1992) to D/H = (2.18 
± °:?2)xlO-5 toward 7^ Vel derived by Sonneborn et al. (2000). 

Hcbrard & Moos (2003) proposed a low value for the present-epoch deuterium abun- 
dance in the Solar neighborhood, based on indirect measurements of D/H from D/0 and the 
constancy of the 0/H ratio (0/H = (3.43 ± 0.15) xlO"^ from Meyer et al. 1998, updated 
/-value). Using more measurements than Hebrard & Moos (2003), Hebrard (2006) proposed 
then D/H = (0.7 ± 0.2)xl0~^ for the present-epoch deuterium abundance. The high D/H 
ratios observed would be due to unknown systematic errors in the determination of A^(HI), 
unknown deuterium enrichment processes and/or local infall of deuterium-rich gas. 

The possibility that deuterium could be depleted in dust grains was first proposed by 
Jura (1982). Draine (2004) and Draine (2006) explored further the theoretical plausibility 
of this process leading Wood et al. (2004) and Linsky et al. (2006) to propose that the 
depletion of deuterium into dust grains leads to different levels of depletion along sightlincs 
with different histories, producing the observed scatter in D/H. In this scenario, the highest 
measurements of the gas-phase abundance of deuterium in the Solar neighborhood place 
a lower limit on the total deuterium abundance, D/H > (2.31 ± 0.24) x 10^^, which one 
can then compare to the primordial value through chemical evolution models. Supporting 
evidence for this model includes correlations of D/H with the depletion of Ti (Prochaska 
et al. 2005), Fe and Si (Linsky et al. 2006), and with the average sightline density (Oliveira 
et al. 2006). 

Galactic evolution models have to reproduce a number of observables such as the metal- 
licity distribution of late type stars, along with predictions of the evolution of deuterium 
from its primordial value to its present day distribution in the Galaxy, and in particular its 
value in the Solar neighborhood. Romano et al. (2006) have addressed the different views on 
the behavior of the D/H ratio by considering what value of D/H is predicted for the Solar 
neighborhood by chemical evolution models that satisfy the majority of the observational 
constraints (such as the G-dwarf metallicity distribution). They find that a modest star 
formation and a continuous infall of unprocessed gas is required to fit all the observational 
data, leading only to a modest decrease of the deuterium abundance from its primordial 
value. According to the models of Romano et al. (2006), which adopt (D/H)p — 2.6x10"^, 
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D/H in the Solar neighborhood should be in the range of (1.4-2.0) xlO~^, implying astration 
factors <1.8. Low D/H ratios, as the one proposed by Hebrard (2006) are ruled out because 
they produced significant disagreements with other observational constraints, in particular 
with the G-dwarf metallicity distribution. 

Testing the two explanations above for the scatter in D/H, requires more measurements 
of D/H and other species that might yield information about the depletion and physical 
conditions along the sightlines, particularly along high column density sightlines. In this 
work we present analyses of the D/H and other ratios along the extended sightlines to 
HD 41161 and HD 53975. Both of these sightlines have total hydrogen column densities 
greater than 10^^ cm^'^, and according to the variable deuterium depletion model such high 
column density sightlines should be biased by cold, not recently shocked gas, where deuterium 
would be depleted in grains. For the HD 53975 sightline we find D/H = (1.02 ± gJ3)xlO-^ 
a value slightly higher but comparable to those found previously for similar sightlines. The 
high D/H ratio along the HD41161 sightline, D/H = (2.14 ± oaI)xW~^, presents the first 
evidence of D/H variations at high -/V(H), indicating that our understanding of the behavior 
of D/H is still incomplete. 

This paper is organized as follows. The targets and their sighthnes are presented in 

§2. §3 describes the observations and data processing, while §4 presents the analyses. The 
abundance of neutral hydrogen along the two sightlines is discussed in §5. The results are 
presented and discussed in §6. Constrains on the /-values of Cll transitions present in the 
spectra of HD 41161 arc derived in §7. Our work is summarized in §8. All uncertainties are 
quoted at the la level unless noted otherwise. 

2. THE TARGETS AND THEIR SIGHTLINES 

The properties of the two stars are listed in Table 1. Below we discuss each star in 
detail. 



2.1. HD 41161 

HD 41161 (08V, Conti & Alschuler 1971) is a field star (Maiz-Apellaniz et al. 2004) 
that probes gas in the direction I = 164.97° and b = +12.89° up to 1253 pc. Jenkins (1978) 
used Copernicus data to derive log iV(0 VI) = 13.29. Using lUE data Shull & van Stccnbcrg 
(1985) derived log iV(Hl) = 20.98 ± 0.07, while Diplas & Savage (1994a) determined log 
A^(HI) = 21.01 ± 0.08 also from lUE data. Taking into account these two determinations 
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of A^(HI) we adopt for the purpose of this study log A^(HI) = 21.00 ± 0.09. 

Dixon et al. (1998) determined the column densities of H2 {J — - 5) along this 
sightline using ORFEUS-1 data. The higher resolution and signal-to-noise of the FUSE data, 
for HD 41161, allow us to determine H2 column densities with uncertainties smaller than the 
~60% (or more) quoted by Dixon et al. (1998). 

2.2. HD 53975 

HD 53975 (07.5V, Conti & Alschuler 1971) is a spectroscopic binary (Garmany et al. 
1980; Gies et al. 1994), member of the CMa OBI association (Humphreys 1978) probing 
gas up to a distance of 1318 pc in the direction I — 225.68° and b — —2.32°. Based on 
observations obtained with Copernicus, Bohlin et al. (1978) determined log A^(HI) = 21.15 
± g:?^. Shull & van Steenberg (1985) used lUE data to derive log A^(HI) = 21.11 ± 0.08, 
while Diplas & Savage (1994a) determined log A^(HI) = 21.10 ± 0.08, also from I UE data. 
We take the weighted average of the independent A^(H I) determinations based on Copernicus 
and lUE (Shull & van Steenberg 1985) data to adopt for this work log A^(HI) = 21.13 ± 
0.06. 



3. OBSERVATIONS AND DATA PROCESSING 

The FUSE observatory consists of four coaligned prime-focus telescopes and Rowland- 
circle spectrographs that produce spectra over the wavelength range 905 - 1187 A, with 
a spectral resolution of ~15 - 20 km s~^ (wavelength dependent) for point sources. Two 
of the optical channels employ SiC coatings, providing reflectivity in the wavelength range 
~ 905 - 1000 A, while the other two have LiF coatings for maximum sensitivity above 
1000 A. Dispersed hght is focused onto two photon-counting microchannel plate detectors. 
With this arrangement of optical channels (LiF 1, LiF 2, SiC 1, and SiC 2) and detector 
segments (lA, IB, 2A, 2B) the FUSE instrument has 8 channels: LiF lA, LiF IB, LiF 2A, 
LiF 2B, SiC lA, SiC IB, SiC 2A, and SiC 2B. Four channels cover the wavelength range 
1000 - 1080 A while two channels each cover the ranges 900 - 1000 A and 1080 - 1180 A. 
Details about the FUSE mission, its planning, and on-orbit performance can be found in 
Moos et al. (2000) and Sahnow et al. (2000). 

Table 2 summarizes the FUSE observations of HD 41161 and HD 53975. HD 41161 was 
observed by F^^-S through the large (LWRS, 30" x 30") and medium-sized (MDRS, 4" x 20") 
apertures. During the short (58 s) LWRS observation data was obtained with the LiF and 
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SiC channels. During the longer (6520 s) MDRS observation only SiC data was obtained 
due to the flux of this star being very close to the FUSE bright object limit (10~^° erg s~^ 
cm~^). The lower reflectivity of the SiC channels compared to that of the LiF channels 
allows some bright targets to be observed, without damaging effects on the detectors. 

The two-dimensional FUSE spectra are reduced using the CalFUSE pipeline v3^. The 
processing includes data screening for low quality or unreliable data, thermal drift correc- 
tion, geometric distortion correction, heliocentric velocity correction, dead time correction, 
wavelength calibration, detection and removal of event bursts, background subtraction, and 
astigmatism correction. The spectra are aligned by cross-correlating the individual exposures 
over a short wavelength range that contains prominent spectral features and then coadded 
by weighting each exposure by its exposure time, using the CORRCAL software developed 
by S. Friedman. 

The S/N of the LWRS data obtained for the HD 41161 sightline is inferior to that of 
the MDRS data. In our analysis we use SiC data from the MDRS observation (no LiF data 
was obtained during this observation) and only LiF data from the LWRS observation. The 
coadded MDRS data has a signal-to-noise ratio per pixel of 16.87 at ~918.5 A in the SiC 
IB channel and of 32.7 at ~1032.2 A in the SiC lA channel. 

HD 53975 was observed by FUSE for 482 s through the MDRS aperture. Similarly to 
HD 41161 only SiC data were obtained during this observation due to the flux of this star 
being close to the FUSE bright object hmit. We coadded the data for this target in the 
manner described above. The coadded data has a S/N per pixel of 10.25 at ~ 918.5 A in 
the SiC IB channel and of 16.87 at ~1032.6 A in the SiC lA channel. 

Figure 1 and 2 present the FUSE spectra of the two targets. 

4. ANALYSIS 

We use apparent optical depth, curve of growth, and profile fitting methods (AOD, 
COG, and PF, respectively) to determine column densities along the sight lines, whenever 
possible (see e.g. Oliveira et al. 2003, 2006, for further discussion of these methods). Details 
about the particular fitting routine used here ( Owens.f) can be found in Hebrard et al. (2002) 
and Lemoine et al. (2002). 

Table 3 lists the atomic data and equivalent widths of the transitions used in this study. 



The CalFUSE pipeline reference guide is available at littp://fuse. pha.jhu.edu/analysis/pipeline_reference. html 
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A, C, and P are used to indicate which transitions are used with the AOD, COG, and PF 

techniques, respectively. No equivalent widths were measured for the transitions flagged only 
with 'P' due to blendings or in the case of D 1 because there are not enough unblended DI 
transitions to construct a curve-of-growth. These transitions were only used with the PF 
technique. We use the compilation by Morton (2003) for the atomic data and that of Abgrall 
et al. (1993a) and Abgrall et al. (1993b) for the molecular data. 

4.1. Atomic species along the line of sight to HD 41161 

We use the AOD technique to place lower limits on the column densities of C 1, C 111, 
Nil, Pll, cm, and Arl along this sight line and to determine the column densities of D I, 
Nl, Fell. With the COC method wc derive column densities for NI, Fell and 6ni = 8.3 ± 
0.2 km s-i and bpeii = 9.2 ± J;^ km s^^ 

The PF technique is used in a similar manner in the analyses of both sightlines. We 
fit a single absorption component to the low ionization atomic species D 1, Nl, 01, and 
Fell, and another to the molecular species. HI is included in the fit as well, in a separate 
absorption component, with the sole purpose of modeling the continuum in the vicinity of 
the D I fines. We use a fine spread function with a FWHM of 0.0641 A to fit the FUSE data, 
constant for all channels and wavelengths, corresponding to a spectral resolution of ~20 km 
s~^ (see Wilfiger et al. 2005, for a determination of the LSF across the FUSE bandpass and 
for a comparison between obtained both with free and fixed LSFs). We derive a Doppler 
parameter 6 ~ 10 km s~^ for the absorption component containing DI, NI, 01, and Fell, in 
agreement with b derived for N 1 and Fe 11 with the COC technique. This Doppler parameter 
should not be viewed as the classic b composed of turbulent and thermal terms, but rather 
the result of line broadening due to multiple unresolved components along the sightline. Fits 
to some of the fines used with the PF method are presented in Figure 3. N(D I) along this 
sightline is constrained mostly by the A916.2 DI transition, which is free from blends with 
other lines, making A"(DI) independent of A"(H2). 

Table 4 presents the column densities of the atomic species obtained with the different 
methods as well as the adopted values. In this work we adopt values that are roughly at the 
midpoint of the values obtained with the different methods and uncertainties that include 
the extremes of the values obtained with the different methods. As a general rule we do not 
adopt 1(7 uncertainties smaller that 0.05 dex (~12%) because we feel that adopting smaller 
uncertainties could be misleading, due to the presence of unknown systematic effects. 

O VI is not clearly detected in our data. Using the A1032 O VI transition we place a 3a 
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upper limit on log A^(OVI) of 13.0, while Jenkins (1978) determined log A^(OVI) = 13.29 

using Copernicus data. Since no uncertainties were quoted by Jenkins (1978) and given 
the S/N of the Copernicus data near the OVI A1032 transition it is likely that the two 
determinations are compatible. 

ClI is also detected along the HD 41161 sightline. We use the A1088 ClI transition to 
determine A^(CII) along this sight line with the AOD method. Besides the A1088 transition, 
many other Cll transitions (for which only theoretically derived oscillator strengths are 
available) are detected. The high S/N ratio of the data allows us to place constraints on the 
/-values of several of these transitions. This is discussed in detail in § 7. 

4.2. H2 and HD along the line of sight to HD 41161 

Absorption by molecular hydrogen from the rotational levels J = through J = 5 
and by HD (J = 0) is present in the spectra of HD 41161. We use the AOD technique to 
determine the column densities of HD (J = 0) and H2 (J = 4 and J = 5). 

We fit a single- 6 COG to the measured equivalent widths of HD and H2 ( J = 1 through 
5) for the HD 41161 sight line. Transitions from the rotational level J = arc too blended 
with those of other species to allow measurements of equivalent widths. We determine 
= 7.2 ± km s-\ 

We use the PF technique to determine the column densities of HD and H2 along this 
sightline, by fitting one absorption component simultaneously, to multiple transitions of HD 
and H2. We derive bu^ ~ 7 km s~^. 

Table 5 presents the adopted molecular column densities along this sighthne (following 
the same reasoning above, used to adopt the atomic column densities for this sightline). Our 
results agree with those derived by Dixon et al. (1998) except for J = 2, for which they 
quote logiV(J = 2) = 17.6 ± 0.2, while we find logiV(J = 2) = 17.97 ± 0.10. We find log 
A^(H2) = 19.98 ± 0.08. Using the column densities for the J = and J = 1 levels we derive 
the excitation temperature Tqi = 84 ± 7 K. This temperature is in good agreement with 
that found by Savage et al. (1977) from Copernicus observations of 61 sightlines with log 
A^(H2) > 18.0 (77 ± 17 K) and that derived by Rachford et al. (2002) from a FUSE survey 
of sighthnes with > 1 mag (68 ± 15 K). 

We derive log A^(H2) = 19.98 ± []:[]^, which combined with iV(Hl) from §2.1 leads to 
the total hydrogen column density along this sightline log A'"(H) = log (iV(HI) -|- 2iV(H2)) 
= 21.08 ± 0.08, and to the molecular fraction /h^ = 2A^(H2)/A^(H) = (16.0 ± f^^) %. 
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4.3. Atomic species along the line of sight to HD 53975 

The AOD technique is used to place lower limits on the column densities of CI, CIII, 
Nil, PII, SIII, and Arl and to determine column densities of DI, NI, O VI, Fell, and CI I 
along the HD 53975 sight line. 

We use the COG method to derive column densities of NI and Fell along the HD 53975 
sightline. We derive 6ni = 9.6 ± 0.4 km and ^Feii ^ 12 km s~^. 

With the PF technique we fit a single absorption component to the low ionization 
atomic species DI, NI, 01, Fell, and Mgll, and another to the molecular species. HI is 
included in the fit as well, in a separate absorption component, with the sole purpose of 
modeling the continuum in the vicinity of the D I lines. A separate component is used to 
model the absorption by O VI. We use a line spread function with a FWHM of 0.0641 A to 
fit the FUSE data, constant for all channels and wavelengths, corresponding to a spectral 
resolution of ~20 km s~^. Figure 4 presents fits to some of the lines used to determine the 
column densities along this sightline. We derive 6 ~ 15 km for the absorption component 
containing DI. Similarly to HD 41161, this Doppler parameter refiects the broadening of the 
lines due to multiple unresolved absorption components along the line of sight. From the 
4 DI transitions used to constrain A^(DI) (see Table 3) only the A 916.9 line suffers from 
major blends with H2 (J = 3) whose column density is not well constrained (sec § 4.4 below). 
Removing this D I line from the fits does not affect A^(D I) adopted. Table 6 presents the 
atomic column densities along this sightline obtained with the different methods as well as 
the adopted values (see § 4.1 for how values and uncertainties are adopted). 



4.4. H2 and HD along the line of sight to HD 53975 

Absorption by HD and by molecular hydrogen from the rotational levels J = through 
J = 5 is detected along this sightline. We use the AOD technique to determine the column 
densities of the H2 J = 4 and J = 5 levels. 

Similarly to HD 41161 we fit a singlc-6 COG to the measured equivalent widths of H2(J 
= through 5) and we use the PF method to determine the column densities of the different 
H2 J levels and HD by fitting one absorption component to H2and HD. We find that the 
column densities of the levels J — through J =3 derived with the COG method are at least 
a factor of ~ 2 larger than the ones derived with the PF technique. An H2 absorption model 
containing the column densities and Doppler parameter derived with the single-6 COG leads 
to H2 absorption profiles that are inconsistent with the data and that clearly overestimate 
the true A'"(H2) along this sighthne. For this reason we adopt for A'^(H2) along this sightline 
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the results derived with the AOD and PF techniques. Transitions from the J — 2 and J = 
3 rotational levels are on the flat part of the curve of growth, hence the associated column 
densities are very uncertain and are flagged with ':' in Table 5. 

We note that A^(D I) is not affected by the adopted A^(H2) as fits performed where the 
DI hues blended with H2 were removed from the fit (particularly A916.9) lead to A'"(DI) 
consistent with the value adopted in Table 6. 

Table 5 presents the adopted molecular column densities for this sightline (see § 4.1 for 
how values and uncertainties are adopted). We find log A^(H2) = 19.18 ± 0.04 by adding 
the column densities of the J = and J — 1 levels, which contain most of the H2 along this 
sightline. Using the column densities for these two levels wc derive the excitation temperature 
Tqi = 93 ± 9 K, in agreement with that found by Savage et al. (1977) but slightly higher 
than that found by Rachford ct al. (2002). The total hydrogen column density along this 
sightline is then iV(H) = 21.14 ± 0.06 with /h^ = (2.2 ± []j)%. 



5. NEUTRAL HYDROGEN ABUNDANCES 

5.1. Ar(HI) ALONG THE HD 41161 LINE OF SIGHT 

To check the consistency of A^(HI) determined from the Lya transition in the IJJE 
data, we use some of the HI Lyman transitions in the FUSE bandpass to determine A^(HI), 
with the profile fitting technique. The following transitions are used to constrain A'"(HI): AA 
1025, 972, 920, 918, 917, and 916. No stellar model is used in conjunction with the FUSE 
data. Wc fit an absorption model with one component of H2, one component containing 
DI, 01, Fell, and Arl, and another component with HI. The damped wings of the Ly/3 
transition constrain iV(HI), while the weaker Lyman lines constrain 6hi- In addition to 
fitting the Lyman lines we include also other regions of the spectra containing absorption 
lines of H2 and other atomic species, to constrain the shapes of absorption lines of these 
species that are unresolved from the H I Lyman transitions. Figure 5 presents the fit to the 
Ly/3 and Ly7 H I transitions along the HD 41161 sight line. The fit yields log iV(H I) = 21.13 
± olol and 6hi ~ 14 km s"^, in reasonable agreement with log A'"(HI) = 21.00 ± 0.09 adopted 
(see §2.1). Varying the degree of the polynomials used to define the continua for the Ly/3 
and Ly7 transitions leads to log A^(HI) in agreement with 21.13 ± qIoI- 
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5.2. Ar(HI) ALONG THE HD 53975 LINE OF SIGHT 

Similarly to HD 41161, we check the consistency of A^(HI) derived from the lUE and 
Copernicus data by fitting some of the Ly H I transitions in the FUSE bandpass. We fit an 
absorption model with one component of H2, a second component with D I, N I, O I, and Fe II, 
and a third component containing HI. The following HI Lyman lines are used to constrain 
A^(HI): AA 916.4, 917.2, 919.3, 920.9, 972, and 1025. Figure 6 presents the resulting fits to 
the HI Ly/3 and Ly7 transitions. We derive log A'"(HI) = 21.19 ± 0.01 (formal statistical 
errors) and 6hi ~ 15 km s~^, in agreement with the adopted value log A^(HI) = 21.13 ± 
0.06. 



6. RESULTS AND DISCUSSION 

Table 7 summarizes several ratios along the two sightlines. Below we discuss each ratio 
in detail. 

6.1. D/H ratios 

In the discussion below the D/H ratios are determined using D/H = [iV(D) + iV(HD)]/[iV(H I) 
+ 2xiV(H2) + iV(HD)]. 

6.1.1. D/H along the HD4II6I sightline 

Figure 7 presents the D/H ratio as a function of log A^(H). Open squares represent 
values from the literature and the three new values from Oliveira et al. (2006) (references for 
all the ratios can be found in this paper) , filled circles the ratios derived in this work. The 
Local Bubble D/H ratio and uncertainties are represented by horizontal solid and dashed 
fines. The vertical dashed fine at log A^(H) = 19.2 represents the edge of the Local Bubble, 
while the vertical dashed line at log A^(H) = 20.7 represents the approximate position of 
where a new low and constant D/H regime was proposed to exist by Wood et al. (2004) and 
Linsky et al. (2006). 

There are five sight fines with log A'"(H) > 20.7, with D/H ratios available in the fiterature 
(HD 198965, HD 191877, JL9, LS 1274, and HD 90087; Hoopes et al. 2003; Wood et al. 2004; 
Hebrard et al. 2005). For these sightlines with 20.7 < log iV(H) < 21.20, the weighted mean 
is D/H = (0.86 ± 0.08) xlO-^ HD 41161 has N{R) similar to these sightlines, but a D/H 
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ratio inconsistent by ~ 3(7. To make D/H along the HD 41161 sightline consistent with the 
weighted ratio of the other 5 sighthncs, one would have to increase iV(H) by a factor of ~ 
2.49, or log iV(HI) from 21.00 to ~ 21.44. To test whether D/H along this sighthne could 
be high due to A^(HI) being underestimated we performed fits where log A^(Hl) was forced 
to 21.44. Figure 8 presents the fits to HI Ljf3 and Ly7 with log A^(HI) forced to 21.44. 
The degree of the polynomial used to fit the continuum of these two transitions, is the same 
as the one used for the fit where we determined A^(HI) = 21.13 ± g gg. However, in this 
case, the continuum is forced to be artificially high to accommodate the higher A^(HI). In 
addition, the fit grossly overestimates the absorption on the blue side of the core of the 
Ly lines as shown in the two bottom panels of Figure 8. To further check that log A'^(H I) 
= 21.44 overestimates the true HI absorption we retrieved the lUE Lya spectrum for the 
HD 41161 sightline from the MAST archive (no further processing was applied). Figure 9 
presents the fit to HI Lya forcing log A^(HI) = 21.44. Similarly to the fits to Ly/? and Ly7 
discussed above, forcing such a high value of A^(HI) on the data forces the continuum to be 
artificially high and overestimates the absorption near the core of the Lya line. 

The A^(HI) adopted for this sightline is based on measurements performed by two 
different sets of authors (Shull & van Steenberg 1982; Diplas & Savage 1994b) who obtain 
similar results, although both used the same lUE dataset. In addition, fits of the Ly/3 and 
Ly7 H I transitions produce A^(H I) in agreement with the value adopted in §2.1, while forcing 
A^(HI) to the value required to bring D/H to agree with the other low D/H measurements 
in the literature (discussed above) produces fits that clearly overestimate the H I absorption. 
We conclude then that the high D/H along this sightline is unhkely due to A'"(HI) being 
underestimated . 

Another possibihty to explain the high D/H along this sighthne is that A^(DI) is over- 
estimated. This is probably unlikely since saturation tends to move column densities in the 
opposite direction, i.e., underestimation rather than overestimation. However other factors, 
such as continuum placement, could produce the desired effect. A D I column density of log 
A^(DI) = 16.00 is required to bring D/H along this sightline down to the weighted mean 
derived by the 5 sightlines discussed above. To test whether this A^(D I) value is consistent 
with the data we consider the A916.18 and A919.1 DI transitions. With the AOD method 
we derive log N{D I) = 16.37 ± 0.04 from the A916.2 transition and our fits yield log iV(DI) 
= 16.43 ± 0.03. A column density of log iV(DI) = 16.0 clearly underestimates the true DI 
absorption, as shown in Figure 10 where we overplot the above column density on the A916.2 
transition. 

To further check our A'"(D I) measurement along this sightline we can consider the A919.1 
D I transition. This transition, not used in our analysis of this sightline, is blended with H2 ( J 
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= 4) and according to our model slightly saturated. Using our H2 fit model we carefully 
deblend the D I and H2 transitions. We use then the AOD technique to derive a lower limit 
to N(D I) using the deblended D I A919.1 transition. This exercise yields log A'"(D I) > 16.13, 
indicating that it is unlikely that A^(D I) is overestimated. 

We conclude then that the high D /H ratio along this sighthne is unlikely due to either 
underestimation or overestimation of the H I or D I column densities, respectively. 

According to the variable deuterium depletion model proposed by Linsky et al. (2006), 
sightlines that probe gas with high column densities should be biased by cold, not recently 
shocked gas, where deuterium would be depleted in grains. The high D/H toward HD 41161 
might indicate that this scenario occurs at A^(H) higher than what was prcvioiisly thought or 
that the clumping of low D/H values in the literature has another explanation. There is also 
the possibility that this sightline is dominated by gas which has been shocked and where the 
typical traces of a shock, such as highly ionized species (not detected along this sightline), 
have been erased. In this scenario the timescale to erase the traces of the shock would have 
to be smaller than that for D to deplete into the grains. The presence of cold H2 (Tqi — 84 
± 7 K) along this sighthne does not necessarily imply that the gas is cold and undisturbed 
since H2 is only a minor constituent of the sightline. Finally, we cannot rule out that local 
production of D I of the type proposed by MuUan & Linsky (1999) or Prodanovic & Fields 
(2003) is responsible for all or part of the deuterium enhancement along this sightline. 

To determine the iron depletion along this sightline, D(Fe) = Log (Fe/H)gas - Log 
(Fe/H)0, we use the solar abundance ratio Log (Fe/H)© = —4.55 ± 0.05 from Asplund et al. 
(2005). We find for this sightline, D(Fe) = —1.55 ± 0.11. In a manner similar to some 
other sightlines with high D/H ratios, the Fe/H ratio along this sightline does not fit the 
trend of high Fe/H versus high D/H, displayed for most of the cases shown in Figure 3 of 
Linsky ct al. (2006). There arc however a few sightlines in Figure 3 of the paper mentioned 
above that seem to follow a trend different than most of the points, in that for the same 
Fe depletion a higher D/H ratio is found. The D/H and Fe/H ratios along the sighthne to 
HD 41161 seem to follow this trend as well. Linsky et al. (2006) discuss some of the reasons 
that could lead to this sort of behavior and argue that these sightlines could either contain a 
smaller percentage of carbon grains (thus fewer sites for deuterium to deplete) or could have 
had weak shocks pass by which would have evaporated deuterium from the grain mantles. 
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6.1.2. D/H along the HD 53975 sightline 

For the HD 53975 sightline we derive D/H = (1.02 ± 020) in agreement with the 
weighted mean D/H = (0.86 ± 0.08)xl0-^ quoted above. Although the D/H ratio along 
the HD 41161 sightline indicates that D/H is variable for log iV(H) > 20.7, the D/H ratio 
along the HD 53975 sightline follows the trend displayed by the other other five sightlines in 
the literature with log A'"(H) > 20.7. Taking into account this ratio, along with the ratios 
for the five sighthnes mentioned above, yields a weighted mean D/H = (0.89 ± 0.07) xlO~^ 
(with xl = 0.4 for 5 degrees of freedom; the average of the upper and lower la uncertanties 
of individual ratios are used to determine weighted means in this work). The iron depletion 
along this sightline, D(Fe) = —1.75 ± 0.09, follows the trend of low Fe/H versus low D/H 
displayed in Figure 3 of Linsky et al. (2006). 

6.2. O/H Ratios 

The O/H ratio along the HD 41161 sightline, O/H = (9.12 ± l^^DxlO'^, is more than 3a 
away from O/H = (3.43 ± 0.15)xl0~^ derived for 13 lines of sight probing gas within 1500 
pc (most within 500 pc), with 20.18 < log A^(HI) < 21.28 by Meyer et al. (1998) (updated 
/-value). For HD 53975, O/H = (5.37 ± \il)xlO-\ is ^1.7a away from the Meyer et al. 
(1998) O/H ratio. The FUSE-ha.sed O/H measurements present a larger scatter than those 
based on if^T data (see for example Andre et al. 2003; Cartledge et al. 2004). Figure 11 
presents a compilation of O/H measurements performed with FUSE (from this work and 
the compilation by Oliveira et al. 2006, asterisks and filled circles, respectively) and HST 
data (from Cartledge et al. 2004, diamonds). For the measurements we consider only 

sightlines that probe gas beyond the Local Bubble, i. e., log A^(H) > 19.2, since previous 
measurements by Oliveira et al. (2005) have shown that O/H = (3.45 ± 0.19) xlO~^ in the 
LB. The assumption that all the HST O/H measurements can be described by a single mean 
leads to the weighted mean O/H = (3.51 ± 0.10)xlO~"' with xl = 1-35 for 51 degrees of 
freedom. Assuming that the FUSE measurements can also be described by a single mean 
leads to the weighted mean O/H = (3.56 ± 0.13) x 10""^ but with xl — 5.4 for 18 degrees of 
freedom. The scatter in the FUSE O/H measurements is more pronounced in the range 20.0 
< Log A'"(H) < 20.9, where there are only a few HST measurements available. Hence, the 
scatter in the FUSE measurements might be due to there really being more scatter for lower 
column density sightlines, even if the O/H values for the LB show little scatter (Oliveira 
et al. 2005). 

The larger (than Meyer et al. 1998) sample of sightlines probed by the work of Cartledge 
et al. (2004) contains some sighthnes with O/H ratios ~6xl0~^. For example, HD 91824 
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{I = 285.70°, b = 0.07°) and HD 91983 {I = 285.88°, b = 0.05°) are both members of the 

Carina OBI association (located at d ~ 2.5 kpc) and have 0/H = (6.91 ± 092) ^10"^ and 
(5.25 ± J;2l)xl0-^ respectively (Cartledge et al. 2004, based on i/^T data). HD 90087 (/ = 
285.16°, b = —2.13°) is also part of the Carina OBI association and is ~100 pc away from 
the two other stars. For this sightline the FUSE derived 0/H = (5.8 ± 1.0) x 10^^ (Hebrard 
et al. 2005) is in good agreement with the 0/H ratios derived for the other members of this 
association. 

In view of the Cartledge et al. (2004) measurements, the 0/H ratio along the HD 53975 
sightline does not seem to abnormally high, however the ratio along the line of sight to 
HD 41161 does. A^(0 1) is determined using the only unsaturated 1 transition in the FUSE 
bandpass, A974, which is heavily blended with H2 J = 2 and J — 5. One might suspect 
that an erroneous /-value would lead to an inaccurate A'"(OI). However, previous studies 
by Hebrard et al. (2005) indicate that A^(0 1) derived from this transition is consistent with 
N{0 1) derived both from weaker and stronger O I transitions. Unknown blends with stellar 
lines are also the unlikely cause of a seemingly high A^(0 1) due to the high Vgim of HD 41161. 
We suspect that unknown systematic effects when deriving A^(0 1) from the heavily blended 
A974 1 transition might be responsible for at least part of the scatter in the FUSE based 
0/H measurements. 

6.3. N/H, O/N, and D/O Ratios 

We derive N/H = (8.32 ± f^DxlO'^ and N/H = (5.07 ± HDxlO'^ for the HD 41161 
and HD 53975 sightlines, respectively. The first ratio is consistent, within the uncertainties, 
with N/H = (7.5 ± 0.4) xlO"^ derived by Meyer et al. (1997). The N/H ratio for the 
HD 53975 sightline is ~1.7o' away from the Meyer et al. (1997). 

We derive O/N = 10.96 ± and O/N = 11.74 ± along the HD 41161 and HD 53975 
sightlines, respectively. The O/N ratios for HD 41161 and HD 53975 are more than 30" and 
2.7(7, respectively, away from O/N = 4.1 ± 0.3 derived using the values of 0/H and N/H 
determined by Meyer et al. (1998) and Meyer et al. (1997). This is likely a consequence of 
the high 0/H for the two sightlines, in addition to the low N/H for the HD 53975 sighthne. 

We derive D/O = (2.29 ± lf^)xlO-^ and D/O = (1.91 ± '^o.ll)^'^^'^ along the HD 41161 
and HD 53975 sightlines, respectively. Both of these ratios are consistent within the uncer- 
tainties with the weighted average D/O = (1.75 ± 0.18)xl0~^ derived by Hebrard et al. 
(2005) for the distant sightlines. The low dispersion of D/O at high N(R) has been used by 
Hebrard & Moos (2003) and Hebrard (2006) as an argument against a high D/H ratio being 
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representative of the interstellar medium in the Solar neighborhood. However, as pointed 
out above, the large scatter in the FUSE denved 0/H ratios seems to indicate that A^(OI) 
might suffer from unknown systematic errors, which would affect the derived D/0 ratios. If 
the high value of A^(Ol) measured along the HD 41161 sightline is indeed too high due to 
systematic effects related to using the A974 transition, then the D/0 ratio along this sightline 
would be have to be revised upwards becoming inconsistent with the Hebrard et al. (2005) 
D/0 ratio for distant sighthnes. This would be consistent with our claim, from D/H, that 
the deuterium abundance towards HD 41161 is higher than that along the other log A'^(H) 
> 20.7 sightlines. 

Finally, we would like to point out that there are other sightlines in which high D/H, low 
D/0, and high 0/H ratios are also observed, similarly to the HD41161 sightline. These are 
Feige 110 (Friedman et al. 2002; Hebrard et al. 2005), PG 0038+199 (Williger et al. 2005), 
and LSE44 (Friedman et al. 2006). Discussions about the possible causes of these ratios are 
given in the papers quoted for each sightline. 

7. /-VALUES OF CHLORINE 

Chlorine has a unique chemistry, in which it can react exothermically with molecular hy- 
drogen. Chlorine is ionized in regions where hydrogen is mostly in the atomic form, whereas 
in regions where there is a significant amount of molecular hydrogen, chlorine is mostly neu- 
tral with a small fraction in the HCl molecular form (Jura 1974; Jura & York 1978). CIl 
traces then the cold molecular components along the sight lines. Determining accurate ClI 
column densities requires accurate oscillator strengths; however only a few transitions have 
experimentally determined /-values. The FUSE bandpass contains numerous transitions of 
CI I, many of which are detected along the HD 41161 sight line. In addition, the strong 
A1347 ClI transition is also detected in the lUE data for this star. Schectman et al. (1993) 
determined the /-values of several ClI transitions using beam-foil spectroscopy. While the 
/-value of ClI A1347 had been previously measured experimentally, AA 1088, 1097 had only 
theoretically determined /-values. The work by Schectman et al. (1993) presents then the 
first and only up to now, experimental measurement of the /-values of the A1088, and A1097 
transitions, and the most accurate measurement of the /-value of A1347. 

Using two ClI transitions which have experimentally determined /-values (AA 1088, 
1347) we place constraints on the /-values of other transitions whose currently adopted the- 
oretical /-values (see Morton 2003, and references therein) are obviously incorrect. For two 
transitions that are optically thin, i. e., on the linear part of the COG, one can derive the 
relationship between their equivalent widths, Wx^^/Wx^ = (/ixA^)/(/2xA|). We can then use 
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the wavelengths and experimentally derived /-values (as well as their associated uncertain- 
ties) for the A1088 and A1347 CI I transitions to derive the predicted ratio ( 1^1347 /W^io88)pred • 
Table 8 presents the atomic data and equivalent widths of the CI I transitions detected along 
the HD 41161 sight line as well as 1 o" upper limits on the equivalent widths of non-detected 
transitions. The fourth column in the table lists the references for the /-values. Using the 
data on this table we derive (Vri347/iyio88)pred = 2.90 ± 0.33, under the assumption that 
for the HD 41161 sightline these two transitions he in the linear part of the COG. Using the 
measured equivalent widths and their uncertainties from Table 8 for these two transitions we 
derive (VFi347/l^io88)meas = 2.35 ± 0.41. The agreement between the predicted and measured 
equivalent width ratios indicates that the linear COG assumption is a valid approximation 
for this sightline. Since ClI traces the cooler components of a cloud, where hydrogen is 
mostly in molecular form, it is possible that there are more than one CI I components, not 
resolved by FUSE. However, as shown by Jenkins (1986), there are cases in which the com- 
bined equivalent widths of a population of Gaussian-like interstellar absorption components 
exhibit a COG which closely mimics that of a single, pure Gaussian distribution in velocity. 

Although it seems that the approach above is valid for this sight line we take a con- 
servative viewpoint by assuming that CI I A1347 could be mildly saturated, but that CI I 
A1088 is not. We use then the equivalent width of this transition as well as its uncertainty 
to predict the /-values of the other weaker CI I transitions quoted in Table 8 using the labo- 
ratory measured /loss and its uncertainty. For the transition x the predicted fx is then — 
/1088 X (Wx/Wioss) X {Xioss/^xT- Blending of ClI lines with stellar lines is not a concern for 
this sightline due to the high Vsini of this star (see Table 1). In addition, careful modehng of 
the H2 and HD absorption along this sight line assures us that the CI I lines studied here are 
free from blends with these molecular species. In particular, the A1088 Cll transition which 
falls near a CO transition (A1087.868) is not affected by a blend with this molecular species 
as log A^(CO) < 13.5 along this sightline. The last column in Table 8 presents the predicted 
/-values (and 1 a uncertainties), using the equation above, which can then be compared to 
the /-values from the literature quoted in the second column of this table. 

Our results, presented in the last column of Table 8, indicate that some Cll transitions 
need to have their /-values revised by large factors. For instance, the currently used /-value 
for A1085 is underestimated by a factor of ~38, while that for A1031 is overestimated by a 
factor of ~7. While the aim of this work is not to determine empirically the Cll /-values in 
the FUSE bandpass, the set of revised /-values we present in Table 8 should provide a better 
match to the observations and can be used in the future for comparison with new theoretical 
calculations. For a comprehensive study of the oscillator strenghts of many CI I transitions, 
including several not covered here, we direct the reader to the work by Sonnentruker et al. 
(2006, in preparation). 
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Taking into account the uncertainties in the /-value of 1088 as well as on the measured 
equivalent width and assuming that this transition in on the linear COG we derive log 
A^(CII) = 13.72 ± 0.04. 

There is a weak absorption feature around 976 A (Wx — 19.5 ± 1.6 A) which we suspect 
to be from chlorine, as we have not been able to identify it with other lines of other atomic 
or molecular species. We use the (J = 2) A975.3507 and HD (J = 0) A975.5821 lines 
to derive the wavelength of this unknown transition by assuming that it traces the same 
absorption component as H2 and HD. The FUSE wavelength solution is expected to be 
accurate to about 6 km s~^ which at these wavelengths translates into ~0.02 A. We derive 
then A = 976.17 A for this unknown transition. 



8. SUMMARY 

We have used FUSE data to determine the abundance of several atomic and molecular 
species along the sightlines to HD 41161 and HD53975. Together with A^(Hl) from the 
literature, we derive D/H, N/H, 0/H, and D/0 ratios for these sightlines. We find that while 
the D/H ratio along the HD 53975 sightline is consistent with five previous measurements at 
similar A''(H), the ratio toward HD 41161 is surprisingly high and presents the first evidence 
of D/H variations for sighthnes with high A'"(H). The pattern of low D/H for sightlines with 
large -/V(H), pointed out by Hebrard & Moos (2003) and emphasized by Wood et al. (2004) 
and Hebrard et al. (2005) is based on a small number of sightlines. Taking into account 
that at this high A^(H) there is a selection bias against high D/H measurements (Oliveira 
et al. 2006) it is not surprising that only 1 out of 7 sightlines present high D/H ratios. 
As analyses of D/H along other high column density sightlines are being performed, one 
might expect that future measurements will also produce high D/H ratios, indicating that 
our understanding of the mechanisms responsible for the observed D/H pattern are still 
incomplete. Finally, the large number of ClI transitions detected in the spectra of HD 41161 
allows to place significant constraints on the /-values of several transitions for which only 
theoretical values are available. 

This work is based on data obtained for the Guaranteed Time Team by the NASA- 
CNES-CSA FUSE mission operated by The Johns Hopkins University. Financial support 
to U. S. participants has been provided in part by NASA contract NAS5-32985 to Johns 
Hopkins University. Finacial support to G. Hebrard has been provided by CNES. Based 
on observations made with the International Ultraviolet Explorer, obtained from the Data 
Archive at the Space Telescope Science Institute, which is operated by the Association of 
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Universities for Research in Astronomy, Inc. under NASA contract NAS5-26555. Support 
for MAST for non-HST data is provided by the NASA Office of Space Science via grant 
NAG5-7584 and by other grants and contracts. The profile fitting procedure, Owens.f, used 
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Table 1. Stellar Properties 



Star 


d 


I 


b 


Sp. Type 


^sin i 






(pc) 


n 


n 




(km S~^) 




HD 41161 


1253 


164.97 


+12.89 


O8.0 V 


300 


0.20 


HD 53975 


1318 


225.68 


-2.32 


07.5 V 


163 


0.22 



Note. — All stellar properties are from Diplas & Savage (1994a) 
except for Vsin i which are from Howarth et al. (1997). Distances are 
based on photometry and are probably accurate within ~25%. 



Table 2. Log of FUSE observations 



Star 


Program ID 


Aperture 


Mode 


Time (s) 


Date 


CalFUSE 


HD 41161 


P1021001 


LWRS 


HIST 


58 


2001 Feb 20 


3.0.7 




P1021002^ 


MDRS 


HIST 


6520 


2003 Sept 25 


3.0.8 


HD 53975 


P3032301^ 


MDRS 


HIST 


482 


2004 Feb 04 


3.0.7 



^Only SiC data obtained during these observations. 
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Table 3. Atomic data and equivalent widths (mA) for the hnes used in the analyses^ 



Species 


Wavelength (A) 


Log f A 


HD 41161 


HD 53975 


DI 


916.1785 


-0.28 


(A, P) 


(A, P) 




916.9311 


-0.18 




(P) 




917.8797 


-0.07 


(P) 






919.1013 


0.04 




(P) 




920.713 


0.17 




(A, P) 


CI 


945.1910 


2.16 


89.59 ± 2.55 (A) 


51.37 ± 4.30 (A) 


cm 


977.0200 


2.87 


222.90 ± 6.22 (A) 


311.13 ± 17.61 (A) 


NI 


951.0792 


-0.79 


67.53 ± 2.42 (C, P) 


49.96 ± 4.01 (C, P) 




951.2948 


-1.66 


18.52 ± 1.95 (A, C, P) 


24.68 ± 3.76 (C, P) 




953.4152 


1.09 


126.05 ± 2.37 (C) 


144.00 ± 3.45 (C) 




953.6549 


1.37 


132.32 ± 2.60 (C) 


147.30 ± 3.41 (C) 




959.4937 


-1.30 


28.89 ± 1.26 (C, P) 


22.55 ± 2.39 (A, C, P) 




960.2014 


-1.95 




(P) 




963.9903 


1.08 




154.30 ± 5.03 (C) 




964.6256 


0.88 




128.36 ± 8.29 (C) 




1134.1653 


1.22 


158.26 ± 4.92 (C) 






1134.4149 


1.51 


162.92 ± 5.15 (C) 






1134.9803 


1.67 


180.55 ± 6.35 (C) 




Nil 


1083.994 


2.10 


170.00 ± 4.94 (A) 


244.55 ± 9.95 (A) 


01 


974.070 


-1.82 


(P) 


(P) 


OVI 


1031.9261 


2.13 




91.16 ± 12.84 (P) 




1037.6167 


1.83 




42.92 ± 7.22 (A, P) 


Mgll 


946.7033 


-0.17 




(P) 




946.7694 


-0.48 




(P) 


PII 


961.0412 


2.53 


74.90 ± 2.23 (A) 


88.23 ± 4.74 (A) 


SIII 


1012.4950 


1.65 




86.38 ± 10.13 (P) 


CI I 


1088.0589 


1.95 


44.55 ± 1.97 (A) 


14.16 ± 3.46 (A) 


cm 


1071.0358 


1.21 


16.16 ±0.96 (A) 




Arl 


1066.6599 


1.86 


114.17 ± 3.03 (A) 


158.12 ± 7.86 (A) 


Aril 


919.7810 


0.91 




(P) 


Fell 


926.8969 


0.72 




19.76 ± 2.06 (C, P) 
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Table 3 — Continued 



Species Wavelength (A) 


Log fX 


HD 41161 


HD 53975 


1055.2617 


0.90 


46.34 ± 1.31 (C, P) 


46.33 ± 2.44 (A, C, P) 


1062.1517 


0.49 


17.59 ± 2.17 (C, P) 




1063.9718 


0.60 


20.41 ± 2.19 (C) 


19.18 ± 3.70 (A, C) 


1081.8748 


1.13 


88.26 ± 2.36 (C) 


61.79 ± 3.59 (C) 


1083.4204 


0.48 


27.75 ± 2.22 (C, P) 




1096.8770 


1.55 


93.16 ± 3.03 (C) 




1112.0480 


0.84 


34.34 ± 3.78 (C) 




1121.9749 


1.36 


83.15 ± 5.79 (C) 




1125.448 


1.26 


84.07 ± 4.35 (C) 




1127.0984 


0.50 


25.86 ± 5.16 (A, C) 




1133.6654 


0.80 


50.57 ± 3.73 (A, C, P) 




1142.3656 


0.68 


45.39 ± 4.79 (C, P) 




1143.2260 


1.31 


80.79 ± 4.64 (C,P) 




1144.938 


2.08 


122.30 ± 5.36 (C) 





^A, C, and P indicate which methods were used with each transition and stand for 
apparent optical depth, curve-of-growth, and profile fitting, respectively. 
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Table 4. Atomic Column Densities (Log) for HD 41161 



Species AOD COG PF Adopted 



HI ... ... ... 21.00 ± 0.09' 

DI 16.37 ± 0.04 (A916. 18) ... 16.43 ± 0.03 16.40 ± 0.05 

NI 17.06 ± 0.04 (A951.29) 17.00 ± 0:0| 16.94 ± 0.01 17.00 ± 0.06 

01 ... ... 18.04 ± 0.05 18.04 ± 0.05 

Fell 14.96 d= 0:°^ (A1127) 15.02 ± 0.06 14.94 ± 0.03 14.98 ± 0.06 

CI I 13.72 ± 0.04 (A1088) 13.72 ± 0.04 

CI > 14.13 (A945) 77. 77. > 14.13 

cm > 14.04 (A977) ... ... > 14.04 

Nil > 14.52 (A1083) ... ... > 14.52 

PII > 13.63 (A961) ... ... > 13.63 

cm > 14.02 (A1071) ... ... > 14.02 

Arl > 14.52 (A1066) ... ... > 14.52 



^See ^2.1. 



Note. — All uncertainties are la. 
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Table 5. Molecular Column Densities (Log) 



Species HD 41161 HD 53975 



HD (J = 0) 



H2 (J = 


0) 


...{J = 


1) 


...{J = 


2) 


...{J = 


3) 


...(J = 


4) 


...(J = 


5) 


7V(H2) 




7V(H) 





14.57 ± 0.08 
19.63 ± 0.06 
19.72 ± 0.06 

17.97 ± 0.10 
17.40 ± 0.10 
14.95 ± 0.08 
14.34 ± 0.05 

19.98 ± 
21.08 ± 0.08 



13.42 ± 0.10 

18.79 ± 0.05 
18.95 ± 0.05 

17.8: 
17.7: 
14.33 ± 0.05 

13.80 ± 0.05 
19.18 ± 0.04 
21.14 ± 0.06 



Note. — All uncertainties are la. 
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Table 6. Atomic Column Densities (Log) for HD 53975 



Species AOD COG PF Adopted 



HI 








21.13 ± 0.06^ 


DI 


16.11 ± 0.03 (A920.7) 




16.20 ± 0.02 


16.15 ± 0.07 


NI 


16.73 ± 0.05 (A959.49) 


16.89 ± O;!? 


16.84 ± 0.05 


16.80 ± 0.08 


01 






17.87 ± 0.08 


17.87 ± 0.08 


OVI 


13.97 ± 0.07 (A1031) 




13.96 ± 0.03 


13.96 ± 0.05 


Fell 


14.80 ± [];[]| (A1063.9) 


14.87 ± 


14.84 ± Hi 


14.84 ± 0.05 


Mgll 






16.10 ± 0.05 


16.10 ± 0.05 


ClI 


13.53 ± 0.10 (A1088) 






13.53 ± 0.10 


CI 


> 13.73 (A945) 






> 13.73 


cm 


> 14.22 (A977) 






> 14.22 


Nil 


> 14.73 (A1083) 






> 14.73 


PII 


> 13.62 (A961) 






> 13.62 


SIII 


> 14.54 (A1012) 






> 14.54 


Arl 


> 14.67 (A1066) 






> 14.67 



*See § 2.2. 



Note. — All uncertainties are la. 
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Table 7. Ratios of Column Densities 



Ratio HD 41161 HD 53975 



(D/H)xlO^ 


2.14 


± 


0.51 
0.43 


1.02 


± 


0.23 
0.20 


(N/H)xl05 


8.32 


± 


2.09 
1.76 


5.07 


± 


1.45 
1.21 


(0/H)xlO^ 


9.12 


± 


2.15 
1.83 


5.37 


± 


1.35 
1.14 


(Fe/H)xl06 


0.79 


± 


0.20 
0.17 


0.50 


± 


0.10 
0.08 


(D/O)xl02 


2.29 


± 


0.40 
0.35 


1.91 


± 


0.51 
0.43 


0/N 


10.96 


± 


2.10 
1.85 


11.74 


: ± 


3.36 
2.80 


(HD/H2)xl06 


3.89 


± 


1.11 
0.93 


1.74 


± 


0.48 
0.39 


/h. (%) 


16.0 


± 


4.0 
3.4 


2.2 


± 


0.4 
0.3 


(nil) (cm"-'^) 





.31 







.31 



Note. — All uncertainties are la and D/H 
= [iV(DI) + 7V(HD)]/[iV(HI) + 2xAr(H2) + 
Ar(HD)]. 
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Table 8. Equivalent Widths and Atomic Data for ClI along the HD 41161 Sighthne 



Wavelength (A) 


/ 


Ref. 


Wx (mA) 


/pred. 




969.9195 


1.73x10"^ 


1 


4.38 ± 0.93 


(1.0 ± 0.2)xl0" 


_2 


984.2865'' 


2.44x10"^ 


1, 4 


13.18 ± 1.95 


(2.9 ± 0.5) xlO" 


9 


1004.6776 


1.58x10"^ 


1, 4 


28.63 ± 1.39 


(6.1 ± 0.7) xlO" 




1022.0478 


1.22x10"^ 


1 


9.35 ± 1.24 


(1 9 ± 3)xl0" 


-2 


1027.3386 


6.58x10-=^ 


1 


12.08 ± 2.97 


(2.5 ± 0.7)xl0- 


-2 


1031.5070 


1.51x10"^ 


1, 5 


10.23 ± 3.42 


(2.1 ± 0.7)xl0- 


-2 


1085.3035 


5.81 xlO^^ 


2, 5 


12.32 ± 2.16 


(2.2 ± 0.5)xl0" 


-2 


1088.0589 


(8.1 ± 0.7)xl0-2 


3, 5 


44.55 ± 1.97 






1094.7686 


1.66x10-2 


2,5 


40.21 ± 3.04 


(7.2 ± 0.9)xl0- 


-2 


1347.2396 


(1.53 ± 0.11)xlO-i 


3, 5 


104.9 ± 17.4 






978.5932 


1.85x10-^ 


1,5 


<4.5 


<1.0xl0-2 




1043.9860 


1.81x10-2 


1,5 


<2.7 


<5.3xl0-3 




1079.8821 


6.98x10-3 


2,5 


<4.1 


<7.5xl0-3 




1099.5230 


1.14x10-2 


2, 5 


<10.0 


<1. 8x10-2 




1101.9361 


1.13x10-2 


2,5 


<5.1 


<9.1xl0-3 





*See § 7 for a discussion of how /pred. is calculated. 
''Possibly blended with the A984.323 transition. 



Note. — For the last six transitions we quote \a upper limits on the equivalent 
widths. 

References. — References 1-3 refer to the original papers. References 4 and 5 are 
compilations of atomic data. (1) Kurucz & Peytremann (1975), (2) Biemont et al. 
(1994), (3) Schectman et al. (1993), (4) Morton (1991), (5) Morton (2003). 
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910 920 930 940 950 960 



1.2 




1020 1040 1060 1080 




1100 1120 1140 ^ 1160 1180 

Wavelength (A) 



Fig. 1. — FUSE spectrum of HD 41161. The FC/5'£' channel used for each panel is indicated 

at the top, in the left. The SiC data displayed is from the MDRS observation, the LiF one 
from the LWRS observation (no LiF data obtained in the MDRS observation; see Table 2). 
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3.0 




950 960 970 980 990 1000 




1010 1020 1030 1040 1050 




1050 1060 1070 1080 1090 1100 

Wavelength (A) 

Fig. 2. — FUSE spectrum of HD 53975. The FUSE channel used for each panel is indicated 
at the top, in the center. No LiF data was obtained for this target. 
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917.4 917. B 91B.2 918. 6 973. B 
10 



974.1 



959.6 



974.4 



974.4 





950.8 951.0 951.2 951.4 1054.4 1054.9 1055.4 

Wavelength (A) 

Fig. 3. — Fits to some of the lines used to determine column densities along the HD 41161 
sightline (convolved with the instrument LSF). Absorption by HI is representd by magenta, 
D 1 by red, N I by dark blue, O I by green, and Fe II by cyan. Absorption by H2 is represented 
by black dashed lines. See § 4.1 for discussion. 
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1.5 
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917.4 950.95 
4 



951.25 



915.55 



1.0 

0.5 
0.0 

918.5 919.0 919.5 920.0 




920.2 



921.3 




959.3 
3 

2 



959.8 



9 9 . 3 







SIB 


v \ t 

TV ' / 
\ \ / 



974.0 



974.5 




974.6 



920.9 921.3 1055.0 1055.2 1055.4 1055.6 

Wavelength (A) 

Fig. 4. — Fit of some of the lines used to determine column densities along the HD 53975 
sightline (convolved with the instrument LSF). Colors are the same as in Figure 3. Absorp- 
tion by Aril is represented by a black dotted line. See § 4.3 for discussion. 
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1020 1022 1024 1026 1028 
Wavelength (A) 




a; 



1023 1024 1025 1026 1027 
Wavelength (A) 




970 971 972 973 974 975 
Wavelength (A) 



0.8 

0.6 

0.4 

0.2 
0.0 











\ i \\ '' ' / \ 




- It V ^ ■ 











971.0 971.5 972.0 972.5 973.0 973.5 
Wavelength (A) 

Fig. 5. — Fit to HI Ly/3 (first and second panels) and Ly7 (third and fourth panels) transi- 
tions along the HD 41161 sightline, yielding log iV(HI) = 21.13 ± o g|, in agreement with the 
value adopted in §2.1. The cores of the Ly/? and Ly7 lines are presented in more detail in 
the second and fourth panels, respectively. H2 is represented by dashed lines, 01, Sill, Dl, 
and Fell by dash-dot-dot. HI is represented by a thick dotted line. The two top panels are 
from SIA MDRS, while the bottom two are from S2B MDRS. The small bump near 1025.5 
A is due to airglow emission, and does not affect our measurements. See §5.1 for discussion. 
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971.0 971.5 972.0 972.5 973.0 973.5 
Wavelength (A) 

Fig. 6. — Same as Fig. 5 but for the HD 53975 sightline, yielding log A^(H I) = 21.19 ± 0.01, 
in agreement with the value adopted in §2.2. The two top panels are from SIB, while the 
bottom two are from SIA. See §5.2 for discussion. 
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I I I I I I I I I I I I I I I I I I I I I I I I I I I I 11 I I I I I I I 




18 19 20 21 

Log N(H) 

Fig. 7. — D/H as a function of log A^(H). Open squares represent ratios from the literature 
and from Oliveira et al. (2006) (where references to all the individual ratios presented can 
be found), filled circles the two new ratios derived in this work. Dashed and solid horizontal 
lines mark the Local Bubble D/H ratio (1.56 ± 0.04) x 10^^ derived from a compilation of 
measurements by Wood et al. (2004). Dashed vertical lines mark log iV(H) = 19.2 corre- 
sponding to the edge of the Local Bubble, and log -/V(H) = 20.7 after which Wood et al. 
(2004) and Linsky et al. (2006) proposed that D/H was low and constant. 
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Fig. 8. — Forced fit to H I Ly/? (first and second panels) and Ly7 (third and fourth panels) 
transitions along the HD 41161 sightline, with log iV(HI) = 21.44, required to bring D/H 
along this sightline to agree with D/H for other 5 sightlines with similar N(H). Dashed lines 
represent H2, dash-dot-dot 01, Sill, DI, and Fell. HI is represented by a thick dotted line. 
The two top panels are from SlA MDRS, while the bottom two are from S2B MDRS. The 
cores of the Ly/3 and Ly7 transitions, presented in more detail in the second and fourth 
panels, show that log A'^(HI) = 21.44 clearly overestimates the true HI absorption. See 
§6.1.1 for discussion. 
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Fig. 9. — Forced fit to the HI Lja transition along the HD 41161 sightline, with log A^(HI) 
= 21.44, required to bring D/H along this sightline to agree with D/H for other 5 sightlines 
with similar N(H). Dashed lines represent NI, Si III, and DI. The global fit, including HI, as 
well as the continuum, are represented by thick solid lines. The core of the Lya transition, 
presented in more detail in the second panel clearly shows that log A'"(HI) = 21.44 clearly 
overestimates the true HI absorption. See §6.1.1 for discussion. 
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Fig. 10. — Log A^(DI) = 16.00, required to lower D/H along the HD 41161 sightline, super- 
imposed on DI A916.2. Absorption is represented by short-dashed lines for HI, long-dashed 
hues for DI, and dotted hues for H2. The overall fit and continuum are represented by 
solid hues. This figure shows that log A^(DI) = 16.00 clearly underestimates the true DI 
absorption. Compare this figure to the corresponding panel in Fig. 3 where we derive log 
N{D1) = 16.43 ± 0.03. See §6.1.1 for discussion. 
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Fig. 11. — 0/H as a function of Log A'^(H). Filed circles correspond to FUSEhetsed measure- 
ments from the literature (see for e.g. Oliveira et al. 2006, for a compilation of these values). 
The two new values derived here are represented by asterisks and open squares represent 
HST based measurements from Cartledge et al. (2004). 



